Boundary-induced phase transitions in traffic flow 
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Boundary-induced phase transitions are one of the surprising phenomena appearing in nonequi- 
librium systems. These transitions have been found in driven systems, especially the asymmetric 
simple exclusion process. However, so far no direct observations of this phenomenon in real systems 
exists. Here we present evidence for the appearance of such a nonequilibrium phase transition in 
traffic flow occurring on highways in the vicinity of on- and off-ramps. Measurements on a German 
motorway close to Cologne show a first-order nonequilibrium phase transition between a free-flow 
phase and a congested phase. It is induced by the interplay of density waves (caused by an on- 
ramp) and a shock wave moving on the motorway. The full phase diagram, including the effect of 
off-ramps, is explored using computer simulations and suggests means to optimize the capacity of a 
traffic network. 
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One-dimensional physical systems with short-ranged 
interactions in thermal equilibrium do not exhibit phase 
transitions. This is no longer true if the action of exter- 
nal forces sets up a steady mass transport and drives the 
system out of equilibrium. Then boundary conditions, 
usually insignificant for an equilibrium system, can in- 
duce nonequilibrium phase transitions. Moreover, such 
phase transitions may occur in a rather wide class of 
driven complex systems, including biological and socio- 
logical mechanisms involving many interacting agents. In 
spite of the importance of this phenomenon and a num- 
ber of theoretical studies 0-^] , it has never been observed 
directly. So far only indirect experimental evidence for a 
boundary-induced phase transition exists in older stud- 
ies of the kinetics of biopolymerization on nucleic acid 
templates ||||. 

In the present work we report the first direct obser- 
vation of a boundary-induced phase transition in traffic 
flow. Analysis of traffic data sets taken on a motorway 
near Cologne exhibits transitions from free-flow traffic to 
congested traffic, caused by a boundary effect, viz. the 
presence of an on-ramp. These transitions are character- 
ized by a discontinuous reduction of the average speed 

Vehicular traffic on a motorway is controlled by a mix- 
ture of bulk and boundary effects caused by on- and off- 
ramps, varying number of lanes, speed limits, leaving 
aside temporary effects of changing weather conditions 
and other non-permanent factors. The fundamental char- 
acteristic of the bulk motion is the stationary flow-density 
diagram, i.e. the fundamental diagram, which incorpo- 
rates the collective effects of individual drivers behavior 
such as trying to maintain an optimal high speed while 
keeping safety precautions. The qualitative shape of the 
flow-density diagram j(p) is largely independent of the 
precise details of the road and hence amenable to numer- 



ical analysis using either stochastic lattice gas models or 
partial differential equations . 

A by now well-established lattice gas model for traffic 
flow, the cellular automaton model by Nagel and Schreck- 
enberg flo|l, reproduces empirical traffic data |ll[] rather 
well. Fig. [j] shows simulation data for the fundamental 
diagram obtained from the Nagel-Schreckenberg (NaSch) 
model. This has to be compared to measurements of the 
flow j taken with the help of detectors on the motorway 
Al near Cologne which show a maximum of about 2000 
vehicles/hour at a density of about p* = 20 vehicles per 
lane and km Q . At densities below p* one observes free 
flow, while for larger densities one observes congested 
traffic. 
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FIG. 1. Fundamental diagram (flow-density-relation) as 
modeled by the NaSch model with v max = 4, p = 0.25, time 
step =1.0 sec, lattice spacing =7.5 m. The system has 3200 
sites, and the flux is averaged over 10 6 lattice updates. The 
broken and full lines indicate the slope which defines the col- 
lective velocity of spontaneous local traffic jams and the shock 
velocity, respectively. 

In addition to the density dependence of the flow two 
important characteristics are derived directly from the 
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fundamental diagram: the shock velocity of a 'domain 
wall' [|l2| between two stationary regions of densities 
P~,P + 



^shock 
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obtained from mass conservation, and the collective ve- 
locity 



dj(p) 
dp 



(2) 



which is the velocity of the center of mass of a local 
perturbation in a homogeneous, stationary background 
of density p. Both velocities are readily observed in 
real traffic. The collective velocity v c describes the up- 
stream movement of a local, compact jam. In the density 
range 25 . . . 130 cars/km v c ranges from approximately 
— 10 km/hour to —20 km/hour (Fig. [I]) which has to be 
compared with the empirically observed value v ~ — 15 
km/hour The shock velocity is the velocity of the 

upstream front of an extended, stable traffic jam. The 
formation of a stable shock is usually a boundary-driven 
process, caused by a 'bottleneck' on a road. Bottlenecks 
on a highway arise from a reduction in the number of 
lanes and from on-ramps where additional cars enter a 
road 

The experimental data considered here (see Fig. || for 
the relevant part of the highway) show boundary effects 
caused by the presence of an on-ramp. Far upstream 
from the on-ramp, free flow of vehicles with density p~ 
and flow j~ = j(p~) is maintained. Just before the on- 
ramp the vehicle density is p + with corresponding flow 
j + = j(p + ). Note that no experimental data are available 
for p~ , j~ and p + , j + as well as the activity of the ramp. 
The only data come from a detector located upstream 
from the on-ramp [ fl6| which measures a traffic density p 
and the corresponding flow j. 
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FIG. 2. Schematic road design of a highway with an 
on-ramp where cars enter the road. The arrows indicate the 
direction of the flow. The detector measures the local bulk 
density p and bulk flow j. 

Next the effects of the on-ramp are considered. Cars 
entering the motorway cause the mainstream of vehicles 
to slow down locally. Therefore, the vehicle density just 
before the on-ramp increases to p + > p~ . Then a shock, 



formed at the on-ramp, will propagate with mean veloc- 
ity v shock (see ([j])). Depending on the sign of v shock , two 
scenarios are possible: 

1) Vshock > (i.e. j + > In this case the shock prop- 
agates (on average) downstream towards the on-ramp. 
Only by fluctuations a brief upstream motion is possi- 
ble. Therefore the detector will measure a traffic density 
p = p~ and flow J = j~ ■ 

2) v s hock < (i.e. j + < The shock wave starts prop- 
agating with the mean velocity v shock upstream, thus ex- 
panding the congested traffic region with density p + . The 
detector will now measure p — p + and flow j = j + . 

Let us now discuss the transition between these two 
scenarios. Suppose one starts with a situation where 
j + > j" is realized. If now the far-upstream-density 
p~ increases it will reach a critical point p cr n < p* above 
which j~ > j + , i.e., the free flow upstream j~ prevails 
over the flow j + which the 'bottleneck', i.e. the on-ramp, 
is able to support. At this point shock wave velocity 
Vshock will change sign (see (0)) and the shock starts trav- 
eling upstream. As a result, the stationary bulk density 
p measured by the detector upstream from the on-ramp 
will change discontinuously from the critical value p cr it 
to p + . This marks a nonequilibrium phase transition of 
first order with respect to order parameter p. The discon- 
tinuous change of p leads also an abrupt reduction of the 
local velocity. Notice that the flow j = j + through the 
on-ramp (then also measured by the detector) will stay 
independent of the free flow upstream from the congested 
region j~ as long as the condition j~ > j + holds. 

Empirically this phenomenon can be seen in the traffic 
data taken from measurements at the detector Dl on the 
motorway Al close to Cologne 0. Fig. ^ shows a typical 
time series of the one-minute velocity averages. One can 
clearly see the sharp drop of the velocity at about 8 a.m. 
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FIG. 3. Time series of the velocity. Each data point repre- 
sents an one-minute average of the speed. Shown are empiri- 
cal data (from [7]) in comparison with the results of computer 
simulations of a simplified model (see text). 

Also the measurements of the flow versus local den- 
sity, i.e. the fundamental diagram (Fig. |J), support our 
interpretation. Two branches can be distinguished. The 
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increasing part corresponds to an almost linear rise of 
the flow with density in the free-flow regime |13|. In ac- 
cordance with our considerations this part of the flow 
diagram is not affected by the presence of the on-ramp 
at all and one measures j = j~ which is the actual up- 
stream flow. The second branch are measurements taken 
during congested traffic hours, the transition period be- 
ing omitted for better statistics. The transition from free 
flow to congested traffic is characterized by a discontinu- 
ous reduction of the local velocity. However, as predicted 
above the flow does not change significantly in the con- 
gested regime. In contrast, in local measurements large 
density fluctuations can be observed. Therefore in this 
regime the density does not take the constant value p + as 
suggested by the argument given above, but varies from 
20 veh/km/lane to 80 veh/km/lane (see Fig. |). 

One should stress here that congested traffic data are 
usually not easy to interpret, because the traffic condi- 
tions (mean inflow and outflow of cars on the on- and off- 
ramps, and so the bulk mean flow) are changing in time. 
According to our arguments, in a congested regime the 
detector measures j + , solely due to the on-ramp activity. 
Therefore, j + (t) = j(p + (t)) < j~(t) must be satisfied. 
During times of very dense traffic one expects always 
cars ready to enter the motorway at the on-ramp, thus 
guaranteeing a sufficient and approximately constant on- 
ramp activity. The measured flow is constant over long 
periods of time which is in agreement with the notion 
that the transition is due a stable traffic jam. Sponta- 
neously emerging and decaying jams would lead to the 
observation of a non-constant flow. 
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FIG. 4. Measurements of the flow versus local density be- 
fore an on-ramp on the motorway Al close to Cologne. The 
detector is located at a distance 1 km upstream from on-ramp. 
Data are collected over a period of 12 days. 

The use of our approach is not limited to a qualitative 
explanation of the traffic data. Beyond that it can also 
be used to calculate the phase diagrams of systems with 
open boundary conditions for a large class of traffic mod- 
els. We modeled a section of a road with on-ramp on the 
left and off-ramp (on-ramp) on the right using the NaSch 



cellular automaton [|TC)|]. We modify the basic model by 
using open boundary conditions with injection of cars 
at the left boundary (corresponding to in-flow into the 
road segment) and removal of cars at the right bound- 
ary (corresponding to outflow). Therefore it can also be 
regarded as a generalization of the asymmetric simple 



exclusion process |17 to particles with higher velocity. 

During the simulations local measurements of the ve- 
locity have been performed analogous to the experimen- 
tal setup. For comparison the results of the computer 
simulations have been included in Fig. 0. Note that even 
the quantitive agreement with the empirical data is very 
good. This has been achieved by using a finer discretiza- 
tion of the model, i.e. the length of the cell is considered 
as I = 2.5m. The results were obtained for L = 960, 
p = 0.25 and v max — 13. We kept the input probability 
a = 0.65 constant. Then the free-flow part is obtained 
using (3—1. and the congested part using f3 — 0.55. The 
transition was observed at ten minutes after we reduced 
the output probability. The "detector" was located at 
the link from site 480 to 481. 

Fig. H shows the full phase diagram of the NaSch model 
with open boundary conditions. It describes the station- 
ary bulk density p as a function of the far-upstream in- 
flow boundary density p~ and the effective right bound- 
ary density p + . For the case of an on-ramp (or shrinking 
road width etc.) at the right boundary corresponds to the 
situation discussed above. Here, the density is increased 
locally to p + > p~ . In agreement with the empirical 
observation we find a line of first order transitions from 
a free flow (FF) phase with bulk density p = p~ to a 
congested (CT) phase with p = p + . On this line v s hock 
changes sign. 

The case of an off-ramp (or expansion of road space 
etc.) leads to a local decrease p + < p~ of the density. 
Here the collective velocity v c (|2]) plays a prominent role. 
As long as v c is positive (i.e. in the free-flow regime 
p~ < p* , see Fig. |l|), perturbations caused by a small 
increase of the upstream boundary density p~ gradually 
spread into the bulk, rendering p = p~ (FF regime). 
At p~ = p* , v c changes sign and an overfeeding ef- 
fect occurs: a perturbation from the upstream boundary 
does not spread into the bulk ^jji) and therefore further 
increase of the upstream boundary density does not in- 
crease the bulk density. The system enters the maximal 
flow (MF) phase with constant bulk density p = p* and 
flow j(p*) = jmax- The transition to the MF phase is of 
second order, because p changes continuously across the 
phase transition point. 

The existence of a maximal flow phase was not empha- 
sized in the context of traffic flow up to now. At the same 
time, it is the most desirable phase, carrying the maxi- 
mal possible throughput of vehicles j max . For practical 
purposes our observations may directly be used in order 
to operate a highway in the optimal regime. E.g. the flow 
near a lane reduction could be increased significantly if 
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FIG. 5. Phase diagram of the NaSch model with open 
boundaries for p — 0.25, v max — 4. Cars enter the road from a 
reservoir of density p~ , inducing the upstream-density p~ dis- 
cussed in the text. At the right boundary cars leave the road 
into a reservoir of density p+, leading to the on-ramp density 
p + . The solid (dashed) curve denotes the theoretical predic- 
tion for the first (second) order transitions lines obtained from 
the numerically determined flow-density relation. The points 
represent phase transition points for the simulated system of 
size 3200. The phases are: free flow (FF), congested traffic 
(CT), maximal flow (MF) phase. 

the traffic state at the entry would allow to attain the 
maximal possible flow of the bottleneck. This could be 
achieved by controlling the density far upstream, e.g. by 
closing temporally an on-ramp, such that the cars still 
enter the bottleneck with high velocity. 

We stress that the stationary phase diagram Fig. || is 
generic in the sense that it is determined solely by the 
macroscopic flow-density relation. The number of lanes 
of the road, the distribution of individual optimal ve- 
locities, speed limits, and other details enter only in so 
far as they determine the exact values characterizing the 
flow-density relation for that particular road. We also 
note that throughout the paper we assumed the external 
conditions to vary slowly, so that the system has enough 
time to readjust to its new stationary state. Experiment- 
ing with different cellular traffic models in a real time 
scale shows that the typical time to reach a stationary 
state in a road segment of about 1.2 km is of the order 
of 3-5 min, which is reasonably small. Close to phase 
transitions lines, however, where the shock velocity van- 
ishes, this time diverges and intrinsically non-stationary 
dynamic phenomena [ ^9|j20[| take the lead. 

In conclusion, we have shown that traffic data collected 
on German motorways provide evidence for a boundary- 
induced nonequilibrium phase transition of first order 
from the free flowing to the congested phase. The fea- 
tures of this phenomenon are readily understood in terms 
of the flow-density diagram. The dynamical mechanism 
leading to this transition is an interplay of shocks and lo- 
cal fluctuations caused by an on-ramp. Full investigation 
of a cellular automaton model for traffic flow reproduces 



this phase transition, but also exhibits a richer phase di- 
agram with an interesting maximal flow phase. These 
results are not only important from the point of view 
of nonequilibrium physics, but also suggest new mecha- 
nisms of traffic control. 
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